Anomalous Relaxation in the XY Gauge Glass 
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To study relaxation dynamics of the two-dimensional XY gauge glass, we integrate directly 
the equations of motion and investigate the energy function. As usual, it decays exponentially at 
high temperatures; at low but non-zero temperatures, it is found to exhibit an algebraic relaxation. 
We compute the relaxation time r as a function of the temperature T and find that the rapid 
increase of r at low temperatures is well described by r ~ {T — Tg)~^ with Tg = 0.22 ± 0.02 
and 6 = 0.76 ± 0.05, which strongly suggests a finite-temperature glass transition. The decay of 
vorticity is also examined and explained in terms of a simple heuristic model, which attributes the 
fast relaxation at high temperatures to annihilation of unpinned vortices. 
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Recently, dynamics of the physical systems which 
have many metastable states has drawn much inter- 
est In general, such systems show very slow re- 
laxation at low temperatures, where energy barriers sep- 
arating metastable states are sufficiently high compared 
with the thermal energy. The XY gauge glass, which 
has attracted much attention in relation to the vortex- 
glass phase of strongly disordered superconductors [Q, 
is a well-known example of the strongly disordered sys- 
tem with a complex energy-landscape. Such a gauge- 
glass model is of particular interest in two dimensions, 
with regard to the possible glass order at finite temper- 
atures. Although it is known that the two-dimensional 
(2D) XY model with random bond angles does not dis- 
play a finite-temperature Kosterlitz-Thouless (KT) tran- 
sition in the strong-disorder regime corresponding to the 
gauge glass |^ , there has been controversy as to the value 
of the glass transition temperature Tg: Numerical cal- 
culations of the current-voltage {IV) characteristics ^] 
have suggested a zero-temperature transition. However, 
the temperature rangeprobed in Ref. |^ is not suffi- 
cient in view of Ref. [Q, where the IV characteristics 
appear to indicate a glass transition at lower temperature 
Tg w 0.15. Further, the analytical result of the vanish- 
ing glass-order parameter at finite temperatures Q does 
not exclude the possibility of an algebraic glass order, 
characterized by the algebraic decay of the glass corre- 
lation function (see Ref. |^ for the definition). This is 
the counterpart of the usual algebraic order (and the 
associated KT transition) in the 2D pure XY model, 
where the Mermin- Wagner theorem pO| ] does not forbid 
the existence of such a quasi- long range order, i.e., the 
algebraic decay of the spin-spin correlation function. Re- 
cently, a similar point has been made on the nature of 



two-dimensional disordered lattice [Q. Therefore, we 
believe that the question on the finite-temperature glass 
transition is not completely settled. Indeed the possibil- 
ity of algebraic glass order at low temperatures has been 
suggested in Ref. |12|, where it has also been shown that 
the numerical calculation of the defect-wall energy |l^] 
does not rule out the existence of an algebraic glass or- 
der at finite temperatures. 

In this work, we study dynamics of the 2D XY gauge 
glass in the absence of an external current and investi- 
gate the possible glass order at finite temperatures. For 
this purpose, we consider an array of resistively shunted 
Josephson junctions with random magnetic bond angles, 
whose equations of motion are derived from the current 
conservation condition at each node. Via numerical inte- 
gration, we compute the energy function and the vortic- 
ity function, and find, as expected, that both quantities 
decay exponentially at high temperatures. In contrast, 
at low temperatures they are found to exhibit algebraic 
relaxation: Such striking difference in the decay behav- 
iors according to the temperature T strongly suggests 
a dynamical glass transition at a finite temperature Tg. 
The behavior of the energy function allows us to estimate 
the relaxation time r, giving the temperature dependence 
T - {T-Tg)-'' withTg = 0.22 ±0.02 and b = 0.76 ±0.05. 
Whereas r is apparently independent of the system size 
at sufficiently high temperature (T 3> Tg), it increases 
strikingly with the system size as T approaches Tg, man- 
ifesting its divergent behavior ai T — Tg in the thermo- 
dynamic limit. We also present a simple heuristic model 
to explain the decay of the vorticity function, and ar- 
gue that unpinned vortices play an important role in the 
relaxation behavior. 

We begin with an L x L square array of resistively 



shunted Josephson junctions with periodic boundary con- 
ditions in both directions. The net current from grain i 
to grain j is written as the sum of the Josephson current, 
the normal current, and the thermal noise current: 



hi = Ic sin( 



(1) 



where (j^i is the phase of the superconducting order pa- 
rameter at grain i, Ic is the critical current of the junc- 
tion, Vij is the potential difference across the junction, 
and R is the shunt resistance. The thermal noise current 
Fjj at temperature T is assumed to satisfy 
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where (• • •) denotes the ensemble average. To describe 
the gauge glass model, the magnetic bond angles Aij's 
are taken to be quenched random variables distributed 
uniformly in the interval (— tt, tt] . The potential is related 
to the phase via the Josephson relation, d{(j)i — (pj ) / dt — 
2eVij/h. This, together with the current conservation at 
each site, allows us to write Eq. (|^) in the form of a set 
of N{= L^) coupled equations: 



dt 



(j)j) + sin(</)j - (f)j - Ay) + jTjij 



0, (3) 



where the summation is over the nearest neighbors of i 
and we have introduced the dimensionless parameter 
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and rescaled time t in units of ?i/ (2ei?/c)- It is of interest 
here to note that Eq. describes the intrinsic dynam- 
ics of the system, based on the current conservation rule. 
This is in contrast to existing studies of dynamical be- 
haviors, where phenomenological Langevin relaxational 
dynamics or Glauber dynamics has been adopted, start- 
ing from the equilibrium Hamiltonian 

To study dynamical behaviors, we integrate directly 
the equations of motion, given by Eq. with L = 8, 
12, 16 and 24. In particular, we use random initial con- 
ditions with the time step At = 0.05, and compute the 
energy function E{t), the auto-correlation function C{t), 
and the vorticity function v{t): 
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where j) the summation over nearest-neighboring 
pairs, (• • •) denotes the average over both the initial 
conditions and disorder configurations. The summation 
'Y^ is taken in the counter-clockwise direction over the 
bonds surrounding plaquette p, and the phase difference 
— Aij is defined modulo 27r. Since the numerical 
results are found to be insensitive to the disorder con- 
figuration, one realization of A^ is used in most cases. 
The auto-correlation function in Eq. (^) measures the 
time correlation of the glass order parameter; the vor- 
ticity function v(t) in Eq. (|^) measures the average chi- 
rality on each plaquette regardless of the sign, and is 
expected to be proportional to the average number den- 
sity of vortices. In the limit t ^ oo, we expect E(t^ and 
vit) should approach the equilibrium values E^'^ and v^'^ ^ 
respectively, at a given temperature T, while C(t) ap- 
proaches XjN regardless of the temperature (see below). 

We first summarize the relaxation behavior of the en- 
ergy function E(t). Figure |^ shows the decay of E(t) 
(in units of hIc/2ekB) of a 16 x 16 array at tempera- 
tures (a) T = 0.1 and (b) T = 0.6. At T = 0.1 it 
is observed that E{t) does not achieve its equilibrium 
value within the time range t < 10^, which corresponds 
to 2 X 10^ time steps. In particular Fig. |] (a) shows 
that for t > 200 E{t) fits well with the algebraic decay 
E{t) — E'^'^ ~ t~°-^°, which is in contrast with the expo- 
nential relaxation E{t)-E'"i ~ exp(-t/19) in Fig. g (b). 
This qualitative change in the relaxation behavior may 
suggest a finite-temperature glass transition with the 
glass transition temperature Tg lower than 0.6. To de- 
termine Tg more precisely, we investigate the normalized 
energy function 



E{t = 0) - E*"! ' 



(8) 



defined to satisfy E{Q) — 1 and E{t oo) = 0. Since a 
system with a complicated free energy landscape in gen- 
eral possesses various scales of the relaxation time [|4|, 
we write the energy function in the form p4|Jl5|| 



dT'P{T')e 



-t/r' 



(9) 



where P{t') describes the distribution of the relaxation 
time. This gives the average relaxation time r as the 
integral of the energy function: 



dT'P{l 



dt'E{t'). 



(10) 



The obtained behavior of r is shown in Fig. || for var- 
ious sizes L = 12, 16, and 24 It is evident that 
at high temperatures (T > 0.5), r is independent of L 
within error bars. At low temperatures, on the other 
hand, r increases with L. This increase becomes more 
pronounced as the temperature is lowered, suggesting di- 
vergent behavior of t in the thermodynamic limit as T 
approaches Tg k, 0.22 from above. Such a rapid increase 
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of the average relaxation time as T is lowered is very well 
described by the form 

r^lT-T,)-" (11) 

with the glass-transition temperature Tg — 0.22 ± 0.02 
and the exponent h — 0.76 ± 0.05, which is represented 
by the solid line in Fig. |^. The dotted line in Fig. |^ is 
the result of the best fit to the form t ~ T~'^, i.e., with 
Tg = 0: It suggests that the behavior of t is incompatible 
with the zero-temperature glass transition. 

This method of estimating the transition tempera- 
ture from the average relaxation time has been used 
in the study of the three-dimensional Ising spin-glass 
model jl^], where both the correlation length in equilib- 
rium and the average relaxation time in dynamics have 
been shown to diverge at the same transition tempera- 
ture. Indeed a very recent study of the correlation length 
^(T) for the 2D XY gauge-glass model showed that i{T) 
diverges at Tg = 0.22 ± 0.01 |]l2|, which is consistent 
with the value obtained dynamically in this work. Here 
the divergent behavior of the correlation length can be 
manifested by characteristic finite-size effects in the nu- 
merical study performed on a system of finite size: Sup- 
pose that the correlation length ^ diverges at Tg like 
^ ~ (T — Tg)~'^ with the appropriate exponent v. In the 
high-temperature regime, the system size L is in general 



with 4)ij{t) = (t)i{t) — <t>j{t)- At t = 0, we have the auto- 
correlation function C(i = 0) = ^ -)- -^N{N - 1) = 1. 
For t oo, we expect that 4>ij{t) and (j)ij(0) become 
independent of each other: 

(cos (t) COS (0)) (cos (t)) (cos (0)}, (13) 

which, combined with (cos (0)) = upon averaging 
over initial configurations, reveals that C{t oo) — 1/N 
at any temperature. The auto-correlation function ob- 
tained numerically at temperatures T = 0.1 and 0.6 is 
plotted in Fig. H, which shows that approach of C(t) to 
the asymptotic value for T < Tg is qualitatively different 
from that for T > Tg. 

We finally study the decay of the vorticity function 
v(t), and propose a heuristic model to understand the na- 
ture of the transition in terms of annihilation of vortices 
and antivortices. The vortices present in the system may 
be classified into three types with regard to their origins: 
quenched vortices induced by the random bond angle 



much larger than ^, and the system is expected not to dis- 
play appreciable size dependence. Near Tg, on the other 
hand, we have ^ ^ T, which leads to strong finite-size 
effects. 

To examine such behavior, we compute the normal- 
ized energy function E{t) at T = 0.6 and at T = 0.25 
and show the results in Fig. ^ and in Fig. |^, respectively. 
At high temperatures (T — 0.6) Fig. || shows that the 
decay of E{t) does not depend on the system size. In 
contrast, at low temperatures (T — 0.25) the relaxation 
becomes slower as the size is increased, as displayed in 
Fig. ||. This size-dependence becomes more distinct as 
the temperature approaches Tg « 0.22 |l^. These dif- 
ferent size-dependent behaviors of E{t) according to the 
temperature are certainly consistent with the divergence 
of the correlation length at Tg « 0.22, which indeed sup- 
ports that both the correlation length in equilibrium and 
the average relaxation time in dynamics diverge at the 
same transition temperature Tg « 0.22. This estimation 
is somewhat higher than the value Tg « 0.15 estimated 
from the IV characteristics [Q . The very slow relaxation 
at low temperatures makes it difficult to reach the sta- 
tionary state in numerical simulations, and presumably 
causes the discrepancy in determining Tg. 

We next investigate the auto-correlation function de- 
fined in Eq. (^), which can be cast in the form: 



(12) 



(type I), thermally generated vortices (type II), and vor- 
tices introduced by the random initial conditions (type 
III) . It is plausible to assume that the numbers of vortices 
and antivortices are equal for each type. Since vortices of 
type III annihilate as time goes on, only vortices of types 
I and II exist in equilibrium; thermal vortices (of type II) 
exist at nonzero temperatures (T > 0) It is further 
assumed that vortices of type III annihilate with antivor- 
tices of the same type or unpinned antivortices of types 
I and II. Since a vortex and an antivortex attract each 
other, and annihilate when they meet, we suppose that 
a vortex-antivortex pair with the size (i.e., the distance 
between the vortex and the antivortex) smaller than / an- 
nihilates after time t. To obtain the transient behavior of 
the vorticity function, we consider a vortex configuration 
consisting of n antivortices of type III and no unpinned 
antivortices of types I and II. The average distance a be- 
tween the randomly positioned antivortices is given by 
a w y^S/TT{n + no), where S is the total area. We now 
consider a single vortex of type III added at a random po- 
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sition. The probability P of the distance between this ad- 
ditional vortex and the nearest antivortex being smaller 
than I may be written as P ~ ttP /a^ « (7r^Z^/S')(n + no), 
which leads to the decay rate for n vortices of type III 



dn 



2;2 



n{n + no) 



(14) 



Note here that no has been assumed to take its equi- 
librium value and the fluctuations of no have been ne- 
glected. Thus Eq. (14) is valid for t greater than the 
"equilibration" time tg of no beyond which no(t) reaches 
its equilibrium value no. It is straightforward to obtain 
the solution of Eq. (0): n(i) = (no/2) [coth(cii+C2) -1], 
where Ci = nQir'^P /2tS and C2 = (1/2) log[l -I- no/n(f)] — 
cii. The decrease of the vortices of type III results in the 
relaxation of the vorticity function. Therefore n is pro- 
portional to v — v'^'^, and the vorticity function is expected 
to exhibit the behavior 



v{t) - v"'^ ~ [coth(cit + C2) - 1], 



(15) 



again for t is larger than tQ. 

We compute numerically the vorticity function, and ex- 
hibit the results in Fig. | for (a) T = 0.1 and (b) T = 0.6. 
At high temperatures (T = 0.6), Fig. ^ (b) shows that 
the numerical data indeed fit very well with Eq. ( |l5| ) 
with ci « 0.035 and C2 ~ 0.087, even for rather small 
t (> 1). On the other hand, at low temperatures, the 
equilibration time to grows large and is expected to be 
larger than the observation time t, thus making Eqs. ( p^ ) 
and (^5|) invalid. In this case the decay of the vorticity 
function is described by an algebraic function rather than 
Eq. (|l|). In Fig. | (a), the vorticity function at T = 0.1 
is shown to be well described by the algebraic function 
v{t) — v'^'^ ^ We have also measured the average 

relaxation time Tvor of the normalized vorticity function, 
defined similarly to Eq. (H), and found that it is also di- 
vergent near Tg, again supporting a finite-temperature 
glass transition. 

In summary, we have investigated dynamics of a two- 
dimensional XY gauge glass via numerical integration of 
equations of motion obtained from current conservation 
conditions. At low temperatures both the energy func- 
tion and the vorticity function have been found to relax 
algebraically, which is in contrast with the exponential 
decay at high temperatures. The auto-correlation func- 
tion at low temperature has also been shown to relax very 
slowly. The relaxation time r has been computed and 
shown to follow t - (T - Tg)-°-'^^ with Tg = 0.22 ± 0.02. 
Such divergence of the relaxation time in turn implies 
the divergence of the correlation length, or slowly (e.g., 
algebraically) decaying spatial correlations ||l^. Finally, 
the relaxation of the vorticity function has been shown 
to be closely related to diffusion of unpinned vortices. 
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FIG. 2. Relaxation time r versus the temperature T — Tg 
in the log-log scale for system sizes L — 12, 16, and 24 with 
Tg = 0.22 ± 0.02. Whereas t is independent of L at high 
temperatures, it strongly depends on L as the temperature 
is lowered. The rapid increase of r fits well with the form 
(T — Tg)^'''^®, which is represented by the solid line. Inset: 
T versus T in the log-log scale. The dotted line is the re- 
sult of the least-square fit to the form T~°, displaying the 
incompatibility with the zero-temperature glass transition. 



FIG. 1. Relaxation of the energy function at temperatures 
(a) r=0.1 and (b) T=0.6. The system size is 16 x 16 and the 
data have been averaged over more than 10^ samples with 
different disorder configurations and initial conditions, and 
error bars correspond to two standard deviations. In (a) the 
solid line represents the algebraic decay E{t) — E"'' ~ f-o.m 
with E'^'' = —0.767, while the solid line in (b) corresponds 
to E{t) - E"" ~ exp(-t/19) with E"" = 0.593. As the tem- 
perature is decreased, the decay behavior changes from the 
exponential to the algebraic. 
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FIG. 3. Relaxation of the normalized energy function E 
at temperature T = 0.6. It is evident that E{t) is indepen- 
dent of the system size. Typical errors axe quite small, of the 
order of 10"^. 
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FIG. 4. Relaxation of the normalized energy function E 
at temperature T = 0.25 for sizes (a) L = 8, (b) L = 12, and 
(c) L = 16. As the system size is increased, the relaxation 
becomes slower. 
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FIG. 5. Relaxation of the auto-correlation function C{t) 
at temperatures T = 0.1 and 0.6 for 16 x 16 array (A'^ = 256). 
Since C{t) ^ 1/7V as t ^ oo, C'{t) - 1/N is plotted as a 
function of time t. 
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FIG. 6. Relaxation of the vorticity function at tempera- 
tures (a) r=0.1 and (b) T=0.6. The system size and other 
conditions are the same as those in Fig. |l| In (a) the solid 
line represents v{t) — 1.678 ~ t~°''^^ , while the solid line in 
(b) corresponds to v{t) - 1.702 ~ coth(0.035f-|-0.087) - 1 [see 
Eq. Note the difference in the vertical scale. 
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